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Abstract of Thesis Presented to the Graduate School 

of the University of Florida in Partial Fulfillment of 

the Requirements for the Degree of Master of Science 
NON-NEWTONIAN TEMPERATURE AND PRESSURE EFFECTS 


OF A POWDER LUBRICANT SLURRY IN A 
ROTATING HYDROSTATIC STEP BEARING 


By 
William Edward Finn 
May 1993 


Chairperson: Dr. J.E. Peterson 
Major Department: Mechanical Engineering 


The purpose of this research was to investigate the 
pressure and temperature effects of graphite powder 
lubricant added to a carrier fluid of ethylene glycol as 
applied in a rotating hydrostatic step bearing. After the 
testing apparatus was designed and modified, the temperature 
and pressure profiles were determined analytically and 
experimentally. 

The rheological behavior of the non-Newtonian lubricant 
was modeled using a Power Law model developed previously. 


This model has been proven to match experimental data and 


Vesa 8 





was therefore the tool used for comparison with research 
data. 

Ethylene glycol was used as the Newtonian lubricant for 
the initial experimentation, which provided a check on the 
test apparatus and a comparison with the non-Newtonian 
graphite slurry. 

The measured data for the temperature increase 
dependency on rotation for both fluids compared favorably to 
the mathematical predictions. A significantly higher 
temperature difference was seen in the non-Newtonian 
lubricant due to the higher shear rates. 

The pressure profile, while not directly dependent on 
rotation as shown in the mathematical model, demonstrated a 
reduction based upon higher rotation speeds. This reduction 
was greater for the non-Newtonian lubricant. Due to the 
higher temperatures seen in the non-Newtonian lubricant and 
temperature dependent Power Law constants, a pressure loss 
is expected. For the Newtonian lubricant, viscosity is 
dependent upon fluid temperature. 

Therefore, the effects of operating speed and 
temperature on a non-Newtonian lubricant should be 


considered, as well as their greater load-carrying capacity. 
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CHAPTER 1 
BACKGROUND AND SIGNIFICANCE OF RESEARCH 


Background of Hydrostatic Bearings 

Hydrostatic bearings were first demonstrated in 1851 by 
L.D. Girard through an invention involving high pressure 
water-fed bearings used for a railway propulsion system and 
patented in 1865 (1, p. 225; 2, p. 435). This externally 
pressurized or hydrostatic system was further demonstrated 
at the 1878 Paris Industrial Exposition when a heavy block 
with four feet was floated on a thin film of lubricant that 
was pumped down each leg (1, 2). The first analysis of a 
thrust bearing was provided by Lord Rayleigh in 1917 using 
the hydrostatic system and solving the load, flowrate and 
metetazonal torque (1; 2, p. 436). 

Current applications include the analysing magnet of 
the Nuclear Structure Facility at the SERC Daresbury 
Laboratory (1, p. 227), the Halle optical telescope (2, p. 
436), the 210-foot-diameter tracking antenna at the NASA-Jet 
Propulsion Laboratory Deep Space Instrumentation Facility 


(3, p. 33-51), and Denver's Mile High Stadium (4, p. 78). 





All applications involve moving extreme weight with minimal 
effort and in some cases, high precision. In 1972, the 
Institution of Mechanical Engineers published proceedings on 
a conference solely devoted to externally pressurized 
bearings (2, p. 437). Since that time, more effort has been 
spent by industry on the study of the lubricants as most of 
the basic mechanical principles of the hydrostatic bearing 

had been solved. 
Significance of the Research 

The addition of solid lubricants to Newtonian fluids 
has been the subject of extensive research. One major 
reason for the addition of solid lubricants is the greater 
Stability of these mixtures when subjected to extreme 
temperatures. The aerospace industry, military, and NASA 
all report using solid lubricants because of their greater 
load and temperature capabilities (5, p. 5; 6, p. 5, 83-86). 
Many commercial applications exist for a Graphite/Glycol 
lubricant mixture including switch gear equipment and high 
temperature oven bearings (5, p. 85-86). 

The U.S. gas turbine industry has started a new program 
in cooperation with NASA and the Department of Defense 


called the Integrated High Performance Turbine Engine 





Technology (IHPTET) program (7, p. 114). The desired double 
propulsion capability from the large speeds in the gas 
turbine shafts result in high temperatures which have 
further emphasized the requirement for nontraditional 
lubricants. The approximate working temperatures of 


different types of lubricants are as follows (8, p. 2): 


Petroleum Lubricants below 100°C 
Synthetic Lubricants 200 tors 57°C 
Powder Lubricants 425 to 675°C 


The powder lubricant is supplied to the system through 
pressurized air or slurries containing the powder (7, p. 
114). These non-Newtonian lubricants are the subject of 


current research towards more accurate behavior prediction. 





CHAPTER 2 
OBJECTIVE OF THE RESEARCH 


This research investigates the pressure and temperature 
affects of graphite powder lubricant added to a carrier 
fluid of ethylene glycol as applied in a rotating 
hydrostatic step bearing. The graphite powder was provided 
by Union Carbide and the ethylene glycol was provided by 
Fisher Scientific. The methodology included laboratory 
testing and mathematical modeling for comparison with the 
experimental results. 

An existing laboratory test rig was extensively 
modified to impart an angular velocity to a hydrostatic step 
bearing and to measure temperature, pressure, flowrate, 
revolutions per minute (RPM), and film thickness. 
Measurements for the ethylene glycol were obtained for use 
as a check on the testing rig and for establishing baseline 
data. 

The goal was to formulate temperature equations 
dependent upon the rotational velocity and compare this 


affect with the pressure distribution across the bearing 





radius. Potential inertial affects from the step bearing's 
rotation and its effects on the graphite particles in the 
ethylene glycol would also be examined in the pressure 
profile. 

An understanding of these fundamental lubrication 
principles in a non-Newtonian fluid is necessary for 
developing a powder lubricant slurry approach in future 


applications. 





CHAPTER 3 
RESEARCH EQUIPMENT AND METHOD 


Laboratory Test Equipment 


The laboratory testing was an important portion of this 
research. Significant time and effort was spent designing, 
fabricating, and redesigning the modifications to the 
existing test rig. The modified hydrostatic test rig 
(Figure 3-1) contains a flat base and a parallel hydrostatic 
step bearing (items 5 and 6, Figure ee. The bearing 
consisted of an upper surface with an outside diameter of 
four inches and an inside diameter of two inches. It was 
connected to the shaft and gear assembly through a ball 
bearing assembly designed and fabricated to maintain a 
parallel surface to the base while rotating. The resultant 
constant film thickness was measured by three LVDT 
displacement sensors (item 4, Figure 3-1) mounted on 
supports which were placed 120 degrees apart on top of the 
hydrostatic step bearing. Since three points define a 


plane, the use of three sensors was sufficient. The voltage 











2) 
S GE > ae 
ff 











| ay 
T omrel Uy) 
al . \ = 





4. Fine odjusting 
thre2d (48 NF) 


: 2 


Fr=7e 
LYOT 


Botion 


6. Boer sng 


Figure 3-1. Modified Hydrostatic Test Rig 





reading of the LVDT displacement sensors was compared to the 
calibration and provided film thickness. The specifications 


for the LVDT sensors are as follows: 


Make Trans-Tek, Inc. 
Model DC-DC Series 240 
02141-0000 
ineweevelts DC (Vv) Ge eo 30 
Linearity % aS 
Freq. Response 3dB down 140 Hz 
Temperature Range S305, Orie) 12 O0gC 
Resolution Infinite 


The film thickness could be adjusted through the fine 
adjusting thread (item 1, Figure 3-1) that connected the 
screw shaft to the frame (item 2, Figure 3-1). The gear 
assembly was belt driven by a controllable, variable speed 
motor (item 3, Figure 3-1). The RPM was measured directly 
through a Pioneer Digital Photo tachometer, model DT-36. 


The motor specifications are as follows: 


Make Dayton Permanent 
Magnet DC 

Model 2M168B 

HP a2 

RPM 25 

Armature Volts 50 DC 

Armature Current 5. 5 amps 


Two thermisters were installed in the base of the 
bearing surface (item 7, Figure 3-1) and measured the 
temperature of the lubricant entering the recess pocket and 


at the perimeter of the hydrostatic step bearing. The 





measured resistance was compared to the thermister's 
calibration and provided the temperature. 

The lubricant was supplied to the bearing surfaces by 
an external pump through an internal orifice. The pump's 


flow rate was adjustable. It's specifications are as 


follows: 
Make Micropump Corp. 
Model 415 
Motor/Pump Coupling Magnetic Drive 
Maximum Flow Rate 630 ml/min 
Rated Maximum Pressure ZO pS y 


The fluid flow was routed through a rotameter flow 
meter. Comparing the reading with the flow meter's 
calibration provided the flow rate. 

The fluid film pressures were measured with dial 
pressure gages at four locations across the bearing surface 
located: 

1) within the recess pocket 

2) 1/4 inch from the inside pocket 

3) 1/2 inch from the inside pocket 

4) 3/4 inch from the inside pocket 

Research Method 
The experimental test rig was used to determine the 
effects of powder lubricants in a carrier fluid subjected to 


rotation in a hydrostatic step bearing. The testing method 
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is described in Chapter 5. First, the various measuring 
devices were calibrated. 
Uncertainty Analysis and Calibration 

The uncertainty analysis for the flowmeter, thermister, 
and LVDT data was calculated using a least squares fit 
technique to determine the deviation width. Since each 
curve fit was linear, the formula (9, p. 158) used to 
determine the standard deviation is written as equation 
(3-1) where the coefficients A and B are determined through 
Paemrinear regression Curve fit. Further, 0, was multiplied 
by 1.7 in order to obtain 10:1 odds or a 90% probability 


that a measurement will fall within 1.7 Standard deviations. 
2 | a 
Oy = 79 201 -A — Bxi) (25515) 


The rotameter flowmeter was calibrated by adjusting the 
flowrate through the pump, marking the flowmeter reading, 
and timing the fill time of the liquid into a beaker of 
known volume. A linear regression or least-squares fit (9, 
p. 154) for a line of the measured data provided the 


calibration equation: 





nab 


y =5.55x — 205.69 fo. 2) 


where y is the flow rate in ml/min and x is the 
flowmeter reading. The calibration curve is shown in figure 
3-2. The flowmeter uncertainty calculated through equation 


(3-1) for both liquids was +11.8 ml/min. 


( 
= 
> 
E 
Lu 
[=~ 
< 
cr 
= 
O 
ad 
iL 


80 90 100 
FLOWMETER READING 


* Data —— Eq. (32) 


Figure 3-2. Flowmeter calibration curve 





The thermisters were calibrated by placing them ina 
water bath, adjusting the temperature, reading the mercury 
thermometer temperature, and recording the thermisters' 
resistance. The uncertainty of the thermometer using half 


the least count was +0.1°F. Thermisters function ona 





a2 
principle that their resistance is lower at higher 


temperatures. This dependency is logarithmic and shown by 


the equation: 


log(#)=a(4-4) +4 (3-3) 


where R, and T, are the baseline resistance in kOhms 
and temperature in Rankine, respectively. The constants A 
and B were determined through calibration for Ena EUG 
thermisters using linear regression. Rearranging the 
results and solving for resistance as a function of 


temperature provided the calibration equations: 


l l 
3083.39( <1, }-.00092 


R, = 11.247) 10 (3-4) 


l l 
31 12.86 t= |-.00027 


Ro = 11.56} 10 (3-5) 


where R,, T, and R,, T, are the resistances and 
temperatures of thermisters one and two, respectively. The 
calibration curves are shown in figures 3-3 and 3-4. Their 


uncertainty calculated through equation (3-1) was + 1 Ohm 


(.001 kOhms). 
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Figure 3-3. 
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Thermister One calibration curve 
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3-4. Thermister Two calibration curve 
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The LVDTs were calibrated by unscrewing the shaft of 
known adjusting thread through equidistant intervals 
beginning at zero when the hydrostatic step bearing was 
touching the flat base. The voltage reading of each LVDT 
was recorded at these intervals. Linear regression of the 
measured data for each LVDT provided the following 


calibration equations: 


LVDT A: yY4=—148.4x,4+ 19.44 (3-6) 
LVDT B: YR_e=—133.8x~R4+ 19.78 (3-7) 
LVDT Cc: Yo=—97.6Xc+19.97 (3-8) 


where y is the LVDT reading in volts and x is the film 
thickness in inches. The calibration curves are shown in 
Figure 3-5. The uncertainty calculated through equation 
(3-1) was +0.04 volts for LVDT A, 40.07 volts for LVDT B, 
and +0.05 volts for LVDT C. 

The uncertainty for the tachometer was +1 RPM as given 
by the manufacturer. It was verified by repeatedly setting 
the maximum RPM value and noting that the value was 


consistently within +1 of 670 RPM. The uncertainty for the 
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Figure 3-5. LVDT calibration curves 


dial pressure gages was +0.5 psi using half the least count. 

The film thickness, h, was calculated as an average of 
the three LVDT readings. In order to determine the overall 
uncertainty of the film thickness during the testing, 
equations (3-6), (3-7), and (3-8) were rearranged to solve 
for film thickness based upon LVDT voltage. The three 


equations were averaged resulting in 


jy = OPC = —.0022y 4 — .0025yz — .0033yc+.161 (3-9) 
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A root-sum-square method (9, p. 73; 10, p. 15; 11, p. 


254) was employed as shown in equation (3-10) where h is 
_— ae ; 
determined from equation (3-9) and = 1s its partial 
t 


derivative with respect to LVDT A, B, or C. The uncertainty 


of LVDT A, B, or C is @, as derived earlier. 


2 
OO, = (2a | (3-10) 


The overall uncertainty of h, the film thickness, was 
+0.0003 inches. This is important because the film 
thickness in the mathematical modeling of the pressure 
profile is raised to the third power resulting in a three 


times greater uncertainty. 





CHAPTER 4 
MATHEMATICAL MODELING 


The mathematical calculations in this chapter were 
formulated in order to model the temperature change and the 
pressure profile in a hydrostatic step bearing while varying 
the film thickness and revolutions per minute (RPM). The 
comparisons of the measured data with the mathematical 


modeling are discussed in Chapter 5. 


f y-axis 


lubricant in (T,, RB) 





Figure 4-1. Schematic of Hydrostatic Step Bearing 
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The equations derived for the temperature and pressure 


relationships of the hydrostatic step bearing shown in 


figure 4-1 involve some basic assumptions (4, pp. 80-88; 8, 


p. 


nO): 


The flow is laminar. The maximum Reynolds number 
using the greatest flow rate and the radius at the 
step is 25.9. A Reynolds number less than 1000 in 
bearing applications is assumed to be laminar flow 
(4, p. 63; 12, p. 343-346). The flow is pressure 
induced. 

The inertial forces resulting from the rotation of 
the bearing are negligible when compared to the 
viscous shear forces and can be neglected. 

The fluids are incompressible. 

The pressure is a function of the radius, 
neglecting any side leakage. 

The temperature is constant across the film 
thickness. 

All the heat generated appears in the fluid and 
corresponds to its temperature rise. 

The pressure at the outside of the bearing is 


assumed to be ambient. 





ue, 
Temperature Equations 

Because of the different properties of Newtonian and 
non-Newtonian fluids, two different mathematical model 
relations were used. 
Ethylene Glycol Relations 

Pure ethylene glycol is a Newtonian fluid and 
therefore, the following relations apply according to the 
geometry of the test rig as shown in the schematic (figure 
4-1). 


The pump loss is given by the relation (4, p. 86), 
Ey = QAP (4-1) 


where Q is the flow rate in cubic inches per second, 
meets Che pressure difference in psi, and E, is energy in 
inch-pounds per second. The friction loss on the bearing 


due to rotation can be solved beginning with Newton's fluid 


friction law, 


Ty = pe (4-2) 





Z0 


Referring to Figure 4-1 and substituting the film 
thickness, h, for the differential width, this is equivalent 


to Newton's equation of viscous friction drag, 
Vv 
— Ss 4- 
f LA; (4=3)) 


Substituting mr* for A and considering a differential 


annulus results, 
at - 
dF = 7; (20rdr) (aa) 


When one of the contacting surfaces is stationary, 
v=@r, where @ is the angular velocity in radians per 
second. Multiplying both sides by r to determine the torque 


on an elemental annulus of radius r and width dr is 





2TUM 
pe rdr (425) 


rae ay = ; 


The total frictional torque (4, p. 88) becomes 





_ 2mm (R 3,._ 2mnol R4-R? 
M i bee. dr = al F (4-6) 











OA 
where wt is the viscosity in reyns, @ is the angular 
velocity in radians per second, h is the film thickness in 
inches, R, is the radius to the step in inches, and R is the 
radius to the outside edge of the bearing in inches. 
The friction loss is given by the equation (4, p. 88) 


_ 2nMN : 
E = aca 


where E, is energy in inch-pounds per second, N is the 
rotation in revolutions per minute (RPM), and M is the 
torque in inch-pounds. Therefore, the total energy loss 


across the bearing from r=0 to r=R can be expressed as 


follows (4, p. 90), 


Ey = Ey+Ey= 12JOyCAT (4-8) 


where 12J is the conversion factor for foot-pounds to 
British thermal units (Btu), Q is the flow rate of the fluid 
in cubic inches per second, y is the weight density of the 
Fluid in pounds per cubic inches, C is the specific heat of 
the fluid in Btu per pounds-°F, and AT is the temperature 


difference across the bearing in °F. The temperature 
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difference can be solved for as a function of the rotation 
Since the flow rate, pressure difference, and film thickness 
are held constant. The pressure difference from r=0 to r=R 
equals the supply pressure since the pressure at the outside 
of the bearing, P(R), is assumed to be ambient. Therefore, 
substituting the values for E, and E, and rearranging to 


solve for AT results 


A mN?y( RRS | 
i Ss 
AT= IAC 21600JOyCh 





(ee) 


Graphite/Ethylene Glycol Relations 

A mixture of ethylene glycol and graphite powder is a 
non-Newtonian fluid and requires different fundamental 
modeling. Currently, there are three accepted rheological 
models (7, pp. 115-116; 8, p. 15). First, the Bingham model 


is stated: 


TNon-N = To + Koy (4-10) 


where T,.,.» 18 the non-Newtonian shear stress in 


Pascals, T, 1s the yield point in Pascals, K, is the Bingham 
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constant in Pascal-seconds, and y is the shear rate in 
seconds’. This model has shown to deviate significantly 
from measurements at low shear rates, however, it is 
accurate at high shear rates where many lubricant films 
operate (7, p. 115). Second, the Hyperbolic model is 


stated: 


1/2 
TNon-N = | (Koy) =) 2K oto | (4-11) 


where the shear stress, yield point, Bingham constant, 
and shear rate variables are the same as the Bingham model. 
This function is similar to the Bingham model at high shear 
rates and has proven to match better with results at low 
Shear rates (7, p. 116). 


Finally, the Power Law model is stated: 
Biv = Oe (ee) 


where T,,,-.» 1S the non-Newtonian shear stress in 
Pascals, K is the Power Law constant in Pascal-seconds”, y 
is the shear rate in seconds", and n is a nondimensional 


Power Law constant. This model has proven useful for 
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describing shear stress-shear rate relations for 
non-Newtonian fluids without a yield point and is useful for 
mec tlow formulations (7, p. 116). 

Since it has been shown that the Power Law model can 
reasonably model a 1:8 ratio of graphite powder to ethylene 
eo 7, Dowtko; 8, pp. 31-32; 13, pp. 86-87) and the 
mathematical expression is uncomplicated to evaluate, it was 
used in this research. 


Table 4-1. Power Law Constants 


Graphite-Ethylene Glycol (1:8) 





The Power Law constants used were taken from Table 4-1 
(7, p. 116) dependent on the temperature of the fluid. 
The pump loss given by equation (4-1) is still 


applicable. The friction loss on the bearing due to 
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rotation can be determined beginning with the Power Law 


model, equation (4-12), rewritten as follows: 
F n 
t=K(#) (4543) 


This is equivalent to the equation of viscous friction 


drag, 
n 
F=KA(2] (4-14) 


Substituting nmr* for A and considering a differential 


annulus results, 
n 
dF = K(z) (2trdr) (4-15) 


Since one of the contacting surfaces is stationary, 
vV=@r, where @ is the angular velocity in radians per 
second. Multiplying both sides by r to determine the torque 


on an elemental annulus of radius r and width dr results, 


rdF =dM = 2nk( 2) (roman (4-16) 
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The total frictional torque becomes 
” OR 
— @ n+2 
M=2nk(® | i (r”"*\dr 


n + n+ 
M =2nk(® | ae | (4-17) 


n+3 


This is equivalent to equation (4-6) for the Newtonian 
fluid where n=1 and K is the viscosity. A conversion factor 
is necessary Since the Power Law constant, K, has units of 
Pascals versus pounds per square inch (psi). The friction 


loss from equation (4-7) expands to the following: 





n 
_ n?KN{ nN RHR 
Bp= 5 &3 n+3 yeeee 


Therefore, the total energy loss from r=0 to r=R can be 
determined according to equation (4-8). The temperature 
difference is solved as a function of the rotation since the 
flow rate, pressure difference, and film thickness are held 
constant. The pressure difference from r=0 to r=R equals 
the supply pressure since the pressure at the outside of the 


bearing, P(R), is assumed to be ambient. Therefore, 
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substituting the values for E, and E, and rearranging to 


Serve for AT results 


Ps m?KN nN \"{ R*3—R2 
AT = Taye + T240200/07C | 30h m3 | — 


Pressure Equations 


The different properties for Newtonian and 
non-Newtonian fluids required different mathematical model 


relations for each as previously discussed in the 





Figure 4-2. Pressure distribution on a unit volume 
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derivations of equations (4-2) and (4-12). The fluid 
flow in the hydrostatic step bearing is pressure induced. 
Ethylene Glycol Relations 

For the geometry of the hydrostatic step bearing 
(Figure 4-1), the pressure distribution on a unit volume is 
shown in Figure 4-2. 

Summing the forces in the horizontal direction and 
rearranging is 


aP _ dt 
dx dy 


(4-20) 
Substituting equation (4-2), it is a simplified form of 


the Navier-Stokes equations (12, p.343) for pressure induced 


Flow between two parallel surfaces, 


ap _ ,,d*v 
dx *"9y? 


(4-21) 
Integrating twice across the y direction, assuming 
dp/dx is a constant across the film thickness, yields 
| ap 


VY) = 374 + Ciy+ Co (4-22) 
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Evaluating the constants of integration, C, and C,, 


with the boundary conditions that v(0)=0 and v(h)=0 results 


Wy) = £207 — yh) (4-23) 


Solving for the flow rate, Q, is shown 


QO=f) v(y)dyb (4-24) 


imtedruataticwmtt 1S stated (4, p. 81; 7, p. 118; 8, p. 
eZ) 


bh3 dp 
O= tre 


one a5) 


where Q is the flow rate in cubic inches per second, b 
is the plate width in inches, h is the film thickness in 
inches, HW is the viscosity in reyns, and dp/dx is the 
pressure gradient. Equation (4-25) is commonly used in 
radial flow situations. Since ethylene glycol is a 
Newtonian fluid, Replacing b by 2nmr, dx by dr, integrating 
and solving for the pressure at any radius, r, becomes (4, 


pecs, Oo, DPal4> 14, Bp. 7) 
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_ 6uQO,_ R 
J ee 0s (4-26) 


where P is the pressure at any point along the radius 
in pounds per square inch (psi), R is the outside radius of 
the bearing, r is any radius between R, and R, Q is the flow 


rate in cubic inches per second, h is the film thickness in 


inches, and w is the viscosity in reyns. 
Graphite/Ethylene Glycol Relations 

Using the same assumptions as the ethylene glycol 
relations with the exception of equation (4-12) for the 
Power Law Model described earlier, it is stated (7, p. 118; 


weep. 16; 14, p. 7) 


i 2nt+l 


d ni ni n 
o=26(4%)'() * (535) 27 





Replacing b by 2nmr, dx by dr, integrating and solving 
for the pressure at any radius, r, becomes (8, p. 16; 14, p. 


8) 


2n+1 
—_ anny 2ntl 4 a) ad a 
P=0Q K em h sees (4-28) 
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The Power Law constants, K and n, are determined from 
table 4-1 for equations (4-22) and (4-23) dependent upon the 
temperature of the fluid. 

For this research, the pressure values for both fluids 
were nondimensionalized by dividing by the inlet supply 
pressure. The flow rate and film thickness were constant so 
that the RPM could be adjusted in order to observe the 
rotational effects on the pressure at the set points along 


the radius. 





CHAPTER 5 
EXPERIMENTAL RESULTS 


Experimental Method 


The initial testing with the laboratory test rig was 
performed with pure ethylene glycol as the lubricant. This 
data provided a check on the instruments and equipment and a 
comparison for studying the effects of adding graphite 
powder to the ethylene glycol. 

The film thickness was set based upon the LVDT 
readings. The flow rate of the lubricant was adjusted so 
that the maximum pressure of the dial gages would not be 

exceeded. The motor was energized and readings were taken of 
the LVDT voltages, flowmeter, pressure gages, and 
thermisters as the rotation was adjusted from 100 to 670 
RPM. When the Graphite/Ethylene Glycol slurry was the 
lubricant, equivalent film thicknesses as the pure ethylene 
glycol were obtained and the above procedure was performed. 
Temperature Results 
The purpose of the thermister data was to determine the 


relation between the temperature difference at the inlet and 


a2 
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the edge of the bearing and the rotation speed (RPM). This 
experimental data were then compared with the mathematical 
predictions as developed in Chapter 4. The data shown in 
Figures 5-1 and 5-2 closely match the predictions within the 
uncertainty. The repeatability was within experimental 


uncertainty. 
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Figure 5-1. Ethylene Glycol Temperature Profile 
Flow rate, h=0.002: 77 ml/min 
h=0.003: 153 ml/min 
h=0.004: 371 ml/min 
Supply Pressure, h=0.002: 40 psi 
h=C.0¢s 25 psi 
h=0.004: 25 psi 


Figure 5-2 clearly shows that the temperature 
difference at higher RPM for the graphite slurry is almost 


double the temperature difference shown in Figure 5-1 at 
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Figure 5-2. Graphite/Ethylene Glycol Temperature 
Profile 
¢ Flow rate, h=0.002: 46 ml/min 
h=0.003: 114 ml/min 
h=0.004: 136 ml/min 
* Supply Pressure, h=0.002: 40 psi 
h=0.003: 30 psi 
h=0.004: 15 psi 


equal RPM and film thickness for the Newtonian fluid, 
ethylene glycol. It is especially prevalent for the 0.002 
inch film thickness. This demonstrates that the presence of 
the graphite powder, thus making the fluid non-Newtonian, 
greatly increases the temperature change. For a rotating 
bearing with microstructures in the lubricant, the increase 
in temperature and subsequent power loss is expected due to 


an increase in the shear stress (15, p. 106). As discussed 
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in Chapter 4, the Power Law constants are dependent upon 
temperature. The values used for the rotation speed up to 
500 RPM corresponded to temperature readings within the 30°C 
range. When the rotation exceeded 500 RPM, the temperature 
readings were higher and fell within the 40°C range. 
Therefore, the Power Law constants were adjusted resulting 
in the discontinuity shown in the theoretical curve in 
Figure 5-2 at h=0.002 inches. The uncertainty for the data 
points was too small to plot. The uncertainty for the 
mathematical model was determined as discussed in Chapter 3 
with equation (3-2). 
Pressure Results 

The purpose of the pressure data was to determine the 
relationship between the pressure profile and the rotation 
speed (RPM). The inertial and temperature effects on the 
pressure are assumed to be negligible when compared to the 
shear stress. It has already been shown that the addition 
of the graphite particles increases the pressure at lower 
flow rates thus, increasing the load capacity of the 
hydrostatic step bearing (8, p. 33; 14, p. 16; 15, p. 106; 
16, pp. 32-33). The effect of rotation on the pressure 


ratio versus the radius ratio is shown in Figures 5-3 
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through 5-8. The measured data compared closely with the 
Mathematical predictions as developed in Chapter 4 at the 


lower rotation speeds within the uncertainty. 


0 H 4 b4 r H : = ‘ x 
Came oomeoo O65 O07 O75 O8 O85 O98 O95 
r/R 


—— Eq. (426) ™ 100RPM * 250 RPM 
x 400RPM & SOORPM 0 600 APM 





Figure 5-3. Ethylene Glycol Pressure Profile 
¢ h=0.002; Flow rate=77 ml/min 

The repeatability was within experimental uncertainty. 
The uncertainty for the data points was too small to plot. 
The uncertainty for the mathematical model was determined as 
discussed in Chapter 3 with equation (3-2). 

The inlet pressure varied slightly during the course of 
the experimentation and therefore, the pressures were 
nondimensionalized by the appropriate supply pressure 


reading resulting in a pressure ratio of 1 at the step of 
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Figure 5-4. Ethylene Glycol Pressure Profile 


* h=0.003; Flow rate=153 ml/min 
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Figure 5-5. Ethylene Glycol Pressure Profile 


¢ h=0.004; Flow rate=371 ml/min 
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Figure 5-6. Graphite/Ethylene Glycol Pressure 
Profile 


* h=0.002; Flow rate=46 ml/min 
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Figure 5-7. Graphite/Ethylene Glycol Pressure 
Profile 


* h=0.003; Flow rate=114 ml/min 
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Figure 5-8. Graphite/Ethylene Glycol Pressure 
Profile 


* h=0.004; Flow rate=136 ml/min 

the bearing for each RPM value. The remaining pressure 
points were nondimensionalized by the same supply pressure. 

The 100 RPM pressure ratios in Figures 5-3 and 5-5 
varied slightly from the mathematical model, but are within 
the uncertainty. As shown in Figure 5-3, the increased 
speed resulted in a slight drop in the interior pressure 
readings whereas the two exit pressures did not vary 
appreciably. The probable reason for this decrease in 
pressure is due to the temperature dependency of the fluid 
viscosity. As noted in the previous results, the 


temperature of the bearing increased due to the raised 
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frictional loss at the higher rotations. As temperature 
increases, the viscosity becomes smaller resulting ina 
smaller pressure ratio. This decrease in pressure is also 
indicated in Figures 5-4 and 5-5. Preliminary temperature 
measurements of the interior of the bearing showed higher 
temperatures than those observed at the inlet or outlet 
points of the fluid flow emphasizing the presence of a 
greater temperature affect on the interior pressure reading 
than the outlet readings. The non-Newtonian pressure 
behavior shown in Figure 5-6 also indicates a drop in the 
pressure ratio. This decrease in the pressure ratio is 
larger than shown for an equal film thickness and RPM of the 
ethylene glycol in Figure 5-3. These results are further 
supported by comparing Figures 5-7 and 5-8 to their ethylene 
glycol counterparts, Figures 5-4 and 5-5. As shown in Table 
4-1, the Power Law constants for the graphite/ethylene 
glycol slurry are also temperature dependent. Therefore, 
the non-Newtonian graphite slurry is also temperature 
dependent. The graphite slurry demonstrated a higher 
pressure drop than the Newtonian fluid, ethylene glycol, 
because of the higher temperatures shown in Figures 5-1 and 


5-2 and discussed previously. 
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The decrease in pressure 1s significant because it 
translates to a loss of load carrying capability of the 
bearing. The higher load carrying capacity of the 
non-Newtonian lubricant is one of the reasons for their use. 
The pressure is still higher than the Newtonian lubricant, 
but it must be considered because the change in temperature 
experienced at extreme rotations will be substantially 
higher than Newtonian lubricants Peeuieine in a higher 


pressure loss. 





CHAPTER 5 
CONCLUSIONS 


The objective of this research was to determine the 
effects of adding a powder lubricant to a Newtonian carrier 
fluid on the lubrication temperature and pressure profile. 
The effects were shown through result comparisons between 
the pure ethylene glycol and the Graphite/Ethylene Glycol 
Sry . 

A detailed uncertainty analysis of the measuring 
devices was undertaken to determine deviations between 
measured and predicted data. 

Results compared favorably with the mathematical 
models. 

The data presented agreed with the theoretical 
predictions based upon the rheological properties obtained 
from previous papers. These agreements give a high level of 
confidence for the use of the equations derived in this 
paper. 

It was noted that the pressure profile of the graphite 


slurry was affected more by the higher rotations than the 
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pure Newtonian fluid. This affect can be expected due to 
the higher shear rate and subsequent temperature increase. 
Since viscosity of the graphite Slurry is temperature 
dependent and the potential applications involve a 
Significant increase to the rotation speed, it is 
recommended that further research should be conducted to 
Study the interior temperature gains and their affect on the 
pressure profile of the bearing. Measuring the temperatures 
of the interior of the bearing in conjunction with the Power 
Law constants' curve fit as functions of temperature could 
result in a complete understanding of the pressure profile 
of a hydrostatic step bearing subjected to a rotation. 

Preliminary testing showed that the temperature inside 
the bearing is higher than those noted at the inlet or 
outlet of the fluid flow. Therefore, any attempt to 
mathematically model the temperature dependency of the 
viscosity or the Power Law constants would have been 
inaccurate. 

Another factor that could be present at significantly 
higher rotation speeds than those tested in this research is 
possible inertial effects of the graphite powder in the 


Slurry mixture. That effect was not apparent in this 
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research since an increase in the pressure profile would be 
expected for such an occurrence. 

The advantages of the non-Newtonian lubricant are their 
increased viscosity and subsequent higher load carrying 
Capacity and their ability to operate at higher 
temperatures. However, when the bearing is operated at high 
rotation speeds, the higher shear rate creates an increased 
temperature which could have a weakening affect on their 
pressure profile. This affect would limit their load 
Carrying capacity, although it still exceeds the capacity of 
the Newtonian lubricant. This thesis has shown the 
fundamental lubrication principles of a graphite slurry 
non-Newtonian fluid that was pressure induced into a 
rotational hydrostatic step bearing. Further investigations 
into higher rotation speeds and interior temperature gain 
are recommended to more fully understand the non-Newtonian 


behavior of the graphite/ethylene glycol lubricant. 
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APPENDIX SAMPLE UNCERTAINTY CALCULATION 


Determining the uncertainty in the film thickness, 
equation (3-9), was accomplished through the use of equation 


(3-10) which is expanded as follows: 


Substituting the values of the partial derivatives of 


equation (3-9) and the uncertainties of each LVDT results: 


cy, = [(-.0022(+.04))? + (-.0025(4.07))? + (-.0033(4.05))? |? 


®, = +0.00026 = +0.0003 


This is the uncertainty stated on page 16 and 


demonstrates the procedure followed for all the uncertainty 


calculations. 
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